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Abstract 

CDF recently measured the quantity j br (b+ H ) ' horn which we determine the ratio ' 

] to be (2.08^^95) x 10~ 3 . In this note, we show that the ratio ^ obtained by dividing the a(B^) by 

the leading order a(b) is consistent with this derived CDF measurement. We calculate the cross section 
a(B^) using the perturbative QCD fragmentation functions of Braaten, Cheung, and Yuan and the 
^ ' corresponding induced gluon fragmentation functions, with the charm-quark mass m c as a parameter. 

' We also estimate the parameter m c from the CDF data and then predict the production rate at Runll 

in 

1. 

CDF |jj and LEP Collaborations ||] recently published their results in search for the final heavy-heavy quark 
bound state - charmed-beauty meson (B c ). The most impressive is the result by CDF, which established 
qq ■ a signal of 4.8cr (from a null hypothesis), using the semi-leptonic decay channels of the B c meson, £>+ — > 



J/ip£ + v, with I = e, /i. CDF measured the ratio 
j*. Tl = ^fj3 H { { B / + g = 0.13218;^ (stat.) ± 0.031 (syst.) i^i (lifetime) , (1) 



where in BR(£?+ — > J /%pi + v) the branching ratios for e and fi are assumed equal, and the last error comes 
from the error in the measurement of the B c lifetime. Based on the following data from Particle Data Book 



BR(5+ - J/ijK+) = (9.9 ± 1.0) x 10~ 4 , = 0.397±8:8Ji , (2) 

BR(B+ -> J/il)l + u) = 2.5 ± 0.5% (3) 



ti : ' (b) 

' and a theoretical calculation 



we are able to deduce this ratio 

a(B+) _ BR(S+ -> J/ipK+) a{B+) 
<r(b) ~ BR(B C + -> J/ipe+v) a(b) 



xu= (2.08iS:§!) x 10" 3 , (4) 



where the error is obtained by adding the relative errors in quadrature. Note that the ratio ^ in Eq. 
(|J) quoted in Particle Data Book represents the fraction of b that hadronizes into a B + meson, which was 
measured at LEP. This fraction is, to a good approximation, independent of px cuts. Thus the ratio ~ 
that we obtained in Eq. (|J) represents the ratio of the cross section of to the cross section of b under the 
same px cut as the CDF measurement 1Z. This ratio has no direct implication that Bj meson is produced 
directly from 6, which is in contrast to B + meson that B + meson is, in general, assumed coming from the 
fragmentation of b. 
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The purpose of this note is to verify that the ratio in Eq. (||) is consistent with a(B+) calculated using 
the perturbation QCD fragmentation functions for b — > £?+ ||, ^, |?], ^ and the corresponding induced gluon 
fragmentation functions, as well as using the leading order (LO) b-quark production. We shall also obtain 
the range of the parameters involved in the fragmentation functions. Once we obtain the parameters we 
can then predict the production rate for the Runll at the Tevatron, where much higher statistics can be 
accumulated. 

The fragmentation approach employed here is different from the full tree-level a 4 calculation || . Compar- 
ison between these two approaches were made in some of the papers in Refs. Basically, the fragmentation 
approach gives a reasonable approximation to the full tree-level calculation as long as pr > 2Mb c - Never- 
theless, one disadvantage of the full calculation is that higher order effects cannot be easily included unless 
the NLO calculation is performed. On the other hand, using the fragmentation approach some important 
higher order effects can be included, namely, the contribution from gluon fragmentation and the contribution 
from higher orbital states below the BD threshold. We shall show that including these contributions we can 
easily account for the ratio in Eq. (0) without employing extreme parameters. 



2. 



In this section, we remind the readers about the importance of the induced gluon fragmentation. The gluon 
fragmentation function for g — > at the initial scale (heavy quark scale) is 0(a s ) smaller than the heavy 
quark fragmentation function for b — > + . Thus, the main source of gluon fragmentation comes from the 
Altarelli-Parisi evolution of the heavy quark fragmentation function. The Altarelli-Parisi evolution equations 
for the fragmentation functions are 



d 



V 

rX dy 



'dy 

y 

1 dy 



Q f dy f dy 

^Q^D g ^ H (z,n) = J —P g ^- b (z/y,fi)D- b ^ H (y,n) + J —P g ^ g (z/y, /i) D g ^ H (y, (i) 



(5) 



(G) 



where H denotes any (be) states, and Pi->j are the usual Altarelli-Parisi splitting functions. The initial scale 
heavy quark and gluon fragmentation functions are H 



D 



b^bc(n 1 So) 



(z,Mo) 



2a s (2m c ) 2 |i? nS (0)| 2 rz(l - zf 



G 



18(1 - 2r)z 



A)^bc(n3Si)( Z !Mo) 



D 



g^B, 



8l7rm3 (1 - fzf 

(21 - 74r + 68r 2 )z 2 - 2f(6 - 19r + 18r 2 )z 3 + 3f 2 (l - 2r + 2r 2 )z 4 

2a s (2m e ) 2 K g (0)| 2 rz(l - zf \ 

TI-Kml (I -fzf I (6 ' 

3(3 -2r + 4r 2 )z 2 - 2f(4 - r + 2r 2 )z 3 + f 2 (3 - 2r + 2r 2 )z 4 
D g^B* {z, M) = for ( u<2(m b + m c ), 



(7) 



(8) 
(9) 
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where r = m c /(mb + m c ), f = 1 — r, //q = m& + 2m c , and i?(0) is the radial wavefunction at the origin. 
They are the initial boundary conditions to the evolution equations in Eqs. (|j)-(^|)- Here we only give the 
S-wave fragmentation functions, which contribute dominantly to B c production, the P-wave fragmentation 
functions can be found in Ref. H. Nevertheless, P-wave fragmentation functions contribute only at 10% 
level to the total B c production. The less determined parameters in the above functions are |-R(0)| and m c . 
The value for i?(0) can be determined in a potential-model calculation |l0). In Ref. |i?is(0)| 2 ranges 

from 1.5 to 1.7 GeV 3 (the extreme value of 3.2 GeV 3 is not used here.) The fixed input parameters of our 
present calculation are tabulated in Table [l], while m c is chosen as a variable parameter in our calculation, 
because the fragmentation function is very sensitive to m c , which appears as m 3 in the denominator: see 
Eqs. (0H§). Overall, we include all n = 1 S-wave and P-wave, and n = 2 S-wave states, which are below 
the BD threshold jl0|], in our calculation. 





n = 1 


n = 2 




4.9 GeV 


4.9 GeV 


Rns(0) 


1.28 GeV 3/2 


0.99 GeV 3/2 


Hi 


10 MeV 




H' & {m) 


1.3 MeV 




cos 6> mix 


0.999 





Table 1: Input parameters to the perturbative QCD fragmentation functions for n = 1 and n = 2. Hi, H' s (m), cos 
are parameters for P-wave states, see Ref. M. 



3. 

We are ready to compute the ratio a(B^) /<r(b) with a(B+) calculated by the fragmentation approach and 
<r(b) by the LO calculation. The "improved" tree-level cross section for is given by 



er(B+) = J dx 1 dx 2 dz f i {x 1 )f j {x 2 ) 



a(ij -> bX, n) D- b ^ B+ (z, fi) + a(ij -> gX, fi) D B + (z, fi) 



, (io) 



where \x is the factorization scale and is chosen to be /i = ht = yPT B + m b- We called this the improved 
cross section because it includes higher order corrections from gluon fragmentation. For b cross section 
we use the LO calculation. When we calculate the ratio of cross sections, the dependence on factorization 
scale, higher-order QCD corrections, parton distribution functions, and mt> are substantially reduced. We 
anticipate the ratio a(B^)/a(b) calculated at tree-level is reasonably accurate without a NLO calculation, 
providing that the present error on B c production is very large. 

We show the ratio a(B+)/a(b) versus PT nlin (b) for m c = 1.2 — 1.7 GeV in Fig. |l|. We note that this 
ratio increases with PT min (^), due to the induced gluon fragmentation contribution. When PT mia (b) increases, 
the scale of the fragmentation function rises and, therefore, the induced gluon fragmentation function also 
increases. If we did not include the induced gluon fragmentation contribution, the ratio <j(B+)/a(b) would 
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Figure 1: The ratio of a(B^)/a(b) versus PT min cut on 6, calculated by fragmentation approach at the Tevatron: 
s/s — 1.8 TeV. A rapidity cut of \y\ < 1 is imposed. The shaded band is the data in Eq. (0), which is derived from 
the CDF data in Eq. ([[]). 

have been a constant, giving rise to a horizontal line coincide with the lower part of the corresponding curve 
in Fig. [l]. Although the gluon fragmentation probability is much smaller than the b fragmentation, the 
production by gluon fragmentation turns out not negligible, because the amplitude squared of the most 
important subprocess gg —> gg is more than an order of magnitude larger than that of gg — ► bb ^| . Figure [l] 
also shows the sensitivity to m c . 

We put the band of a(B+)/a(b) given by Eq. (g) onto Fi g- 0- We note tnat tne CDF data in Ec I- (0) 
is for i?+ and B + with pr > 6.0 GeV and \y\ < 1 We have to convert this pr requirement on B + and 
.B+ to Pt requirement on 6, because the fragmentation spectrum of b is not monochromatic. The average 
momentum fraction (z) for fragmentation of b into B + and _B+ is about 0.7 — 0.8 at the scale [i w 8 — 10 
GeV. Hence, the pr requirement on b becomes 8 — 9 GeV. From Fig. |l| at around PT mln (b) =8 — 9 GeV, 
the shaded CDF band gives 

m c ~ 1.3 — 1.7 GeV , (11) 



with the central value at about 1.45 GeV. Since the error of the ratio in Eq. (W) is large, the range of m c 
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Figure 2: The ratio of a(B^)/a(b) versus PT min cut on b calculated by fragmentation approach at Run II: ^fs = 2 
TeV. The shaded region corresponds to m c ~ 1.3 — 1.7 GeV with the solid line at m c — 1.45 GeV. A rapidity cut of 
\y\ < 1 is imposed. The factorization scale fi — /it = \Jpt- + m f> f° r ^ ne s °lid line, /i = fir/2 for the dashed, and 
/x = 2/ir for the dot-dashed. 

obtained in Eq. (0) is also very wide. 



4. 

Run II at the Tevatron will be at y/s = 2 TeV with a nominal accumulated luminosity of 2 fb _1 . The 
prediction of a(B+) / a(b) for the range of m c ~ 1.3 — 1.7 GeV obtained above in Eq. ( JTl| ) is given in Fig. || 
(shaded region) with the solid line for rn c = 1.45 GeV. It appears that the ratio predicted at yfs = 2 TeV is 
about the same as at yfs = 1.8 TeV. 

Finally, we also demonstrate the dependence of the ratio on the factorization scale, which appears in the 
running a s , the parton distribution functions, and the fragmentation functions. In Fig. 0, the solid line is 
for the original choice of fj, = fiT = \jpr- b + m 1 > while the dashed is for = /j.t/2 and the dot-dashed for 
/i = 2fiT- The dependence on the scale is smaller than the dependence on m c . 
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To summarize we have obtained the ratio a(B+)/a(b) = ^2.08jrJ;g|j x 10~ 3 from the CDF data in Eq. 
(|l|). We have also verified that the prediction by the perturbative QCD fragmentation approach is consistent 
with the CDF data, with m c ~ 1.3 — 1.7 GeV and the central value at 1.45 GeV. The prediction of the ratio 
at Run II is very similar to that at Run I. 
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